Generation-recombination processes via acoustic phonons in a disorded graphene 
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Generation-recombination interband transitions via acoustic phonons are allowed in a disordered 
graphene because of violation of the energy-momentum conservation requirements. The generation- 
recombination processes are analyzed for the case of scattering by a short-range disorder and the 
deformation interaction of carriers with in-plane acoustic modes. The generation-recombination 
rates were calculated for the cases of intrinsic and heavily-doped graphene at room temperature. 
The transient evolution of nonequilibrium carriers is described by the exponential fit dependent on 
doping conditions and disorder level. The characteristic relaxation times are estimated to be about 
150 - 400 ns for sample with the maximal sheet resistance ~5 kf2. This rate is comparable with the 
generation-recombination processes induced by the thermal radiation. 

PACS numbers: 72.80.Vp, 73.61.-b, 78.60.-b 



I. INTRODUCTION 

Transport [l| and optical [2| properties ofgraphene 
as well as noise phenomena in this material |3| are not 
completely understood for the regime of nonlinear re- 
sponse. The treatment of nonequilibrium carriers re- 
quires not only verification the momentum and energy 
relaxation processes but also understanding of the in- 
terband generation-recombination processes which deter- 
mine electron and hole concentrations far from equilib- 
rium (similar transport conditions take place for the bulk 
gapless materials, see review J4| and references therein). 
Effective interband transitions via optical phonons of en- 
ergy huo take place for the energy of carriers greater than 
hujo/2, see Refs. 5 and 6 where the cases of optical ex- 
citation and heating by dc current were analyzed. At 
lower energies, the generation-recombination processes 
become ineffective because the Auger transitions are for- 
bidden due to the symmetry of electron-hole states [7| 
(c.f. with Q). Since the carrier's velocity v ~ 10 8 cm/s 
exceeds significantly the sound velocity s, the interband 
transitions via acoustic phonons are also forbidden due 
to the momentum-energy conservation laws. Only slow 
generation-recombination processes induced by the ther- 
mal radiation are allowed in a perfect graphene. |9j To 
the best of our knowledge, consideration of a disorder ef- 
fect on the interband transitions via acoustic phonons in 
the low-energy region, e < Tiloq/2, is not performed yet. 
Thus, the evaluation of the generation-recombinaton rate 
caused by the interaction of carriers with the acoustic 
phonon thermostat under violation of the momentum- 
energy laws in a disordered graphene (allowed electron- 
hole transitions are depicted in Fig. 1) is timely now. 

In this paper, the calculations are performed for the 
model of short-range disorder whose parameters are 
taken from the mobility data, [l|, [l(| for samples with 
the maximal resistance of 2 - 6 kO per square. The prob- 
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FIG. 1: Interband generation-recombination transitions via 
acoustic phonons with energies ~ Tiuiac (thick arrows) between 
broadened electron-hole (e - h, shown by grey) states in the 
low-energy region, e < huo/2. Thin lines show the ideal dis- 
persion law. 



ability of electron-hole transitions is expressed through 
the averaged spectral density functions and is calculated 
taking into account the contribution of interband inter- 
ference. Due to slowness of the interband transitions, 
the quasiequilibrium distributions of electrons and holes 
with the same temperature are used for the description 
of a temporal evolution of nonequilibrium concentrations 
of carriers. The electron and hole concentrations are also 
connected through the electroneutrality condition with 
the surface charge controlled by a gate voltage. 

The results were obtained for the cases of intrinsic and 
heavily-doped graphene at temperature T and can be 
briefly summarized as follows. The concentration bal- 
ance equation is written through the chemical potential 
normalized to T and the characteristic rate, which is pro- 
portional to a carrier-phonon coupling and increases with 
temperature as T 2 . The transient evolution of nonequili- 
bruim population can be fitted by an exponential decay 
with the relaxation time 150 - 400 ns at room temper- 
ature and a typical disorder level corresponding to the 



maximal sheet resistance ~5 kO. This time scale appears 
to be comparable to the recombination rate via thermal 
radiation and the mechanism under consideration can be 
verified by temperature and temporal measurements. 

The paper is organized as follows. In the next sec- 
tion we present the basic equations which describe the 
generation-recombination processes under consideration. 
In Sec. Ill we evaluate the generation-recombination 
rates and analyze their dependencies on temperature, 
disorder level, and doping conditions. The last section 
includes the discussion of the approximations used and 
conclusions. In Appendix we consider the generation- 
recombination mechanism caused by the interaction with 
the thermal radiation. 



II. BASIC EQUATIONS 

Temporal evolution of carriers in a random potential, 
which are weakly interacting with the acoustic phonon 
modes, is described by the distribution f at over the states 
\a) with energies e a . An exact with respect to a disordere 
effect kinetic equation takes the form |llj 



dt 



a' 



[W a 'afa't (1 - f at ) 

-W aa ,f at {l-f a >t)]- 



(1) 



The transition probability W aa < is written within the 
Born approximation with respect to the carrier-phonon 
interaction with gth phonon mode of frequency cu q : 
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Here the operator x q determines the carrier-phonon in- 
teraction Y] I Xqbq + H.c. ) where b q is the annihilation 

operator of qth mode and N q is the Planck distribution 
of phonons at the equilibrium temperature T. Note, that 
the transition probabilities W a ' a and W aa > are connected 
by Wa'a = cxp [- (e a - e a >) /T] W aa > and 
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due to the particle conservation law. 

The concentrations of electrons and holes, rtt and fit, 
which are averaged over random disorder (such averaging 
is denoted as (■••)), are given by 
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where L 2 is the normalization area and the step func- 
tion 0(±e) appears due to the symmetry of electron-hole 
spectrum [see Eq. (10) below]. Since effective intra- 
band scattering is caused by the phonon thermostat and 



carrier-carrier interaction, the quasiequilibrium distribu- 
tions over the conduction and valence bands are imposed 
during a short-time scales and below we use 



let 



exp 



exp 
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Due to effective energy relaxation the same temperatures 
are established in both bands. At the same time the elec- 
ton and hole concentrations are determined through the 
different chemical potentials, /i^ and /if, respectively. 
The chemical potentials are connected by the electroneu- 
trality condition n t —lit = n s , where the surface charge 
en s is controlled by the gate voltage, V g , according to 
n s = aV g with a ~ 7.2 x 10 10 cm~ 2 /V written for the 
SiC>2 substrate of thickness 0.3 /im. 

The concentration of electrons is governed by the bal- 
ance equation dn t /dt = (dn/dt) ac with the generation- 
recombination rate 



exp 
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We take into account that the intraband transitions 
(when e, e' > 0) vanish in (6) and transform the tran- 
sition probability as follows 



W(e, e') = — 1^6(6 -e a )S (e' - e a >) W aa , \ . (7) 



Further, we introduce the exact spectral density function 
A E (Ik, l>x') ="£6(e- ea) ^g*. ( 8 ) 

a 

which is determined through the double-row wave func- 
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tion $j x with I =1,2. The column ^x is a solution of 
the eigenvalue problem (h + V x )^!x = £q^x written 
through the single-particle Hamiltonian h and a random 
potential V x . Using the definition (2) one obtains the 
probability W(e,e') as follows 

W(e,e') = -j^J2\ C i\ 2 ( N * + 1 ) 5 ( £ ~ £ '- h ^ ( 9 ) 
q 

x f dx f dx'e lq ( x - x ')tr (A e , (x, x') A e (x', x) 

where q is the in-plane wave vector and |CJ 2 is the 
matrix element of deformation interaction. [10] As a 
result, (dn/dt) ac is expressed through the two-particle 
correlation function. Since the main contributions to 
(6) appears from e =/= e' , this correlation function 

can be decoupled according to lA e i (x, x') A s (x',x) 
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FIG. 2: (Color online) (a) Density of states p e versus energy 
at g =0 (1), 0.15 (2), 0.3 (3) and 0.45 (4). (b) Ratio V s = 

y/7>Jpe versus g for e =20 meV (1), 30 meV (2), 40 meV (3), 
and 50 meV (4). (c) Equilibrium concentration of non-doped 
graphene n eq versus g normalized to tit — 0.52(T/hv) 2 . 



4',Ax4,-Ax, where Ar,Ax = (A £ (x,x')\ is the aver- 
aged spectral function given by 2 x 2 matrix. 

Below, we calculate the probability (9) using the model 
of the short-range disorder described by the Gaussian 
correlator (V^V^) = ^ 2 exp[-(x - x') 2 /2l 2 }, where V is 
the averaged amplitude, l c is the correlation length, and 
the cut-off energy E c = vfi/l c exceeds the energy scale 
under consideration. According to Refs. 12 and 13 the 
retarded Green's function in the momentum representa- 
tion takes form: 
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where P^' - [1 ± (<r • p)/p] /2 are the projection op- 
erators on the conduction (+) and valence (— ) bands, 
<x is the isospin Pauli matrix, and g = (V 2 l c /hv) 2 Tr/2 
is the coupling constant. Here we restrict ourselves by 
the Born approximation when the self-energy contribu- 
tion e(A e + ig) is written through the logarithmically- 
divergent real correction and the damping factor. Note, 
that these corrections vanish at e — > 0. From a compar- 
ision with the mobility data [l|, [l(| one obtains that the 
parameters g ~0.45, 0.3, and 0.15 correspond to the sheet 
resistances ~6, ~4, and ^2 kfl per square, respectively. 
The density of states, p £ — — 4Im^ trGf p is shown in 
Fig. 2a and p £ = p- £ , i. e. the electron- hole symme- 
try is not violate due to disorder. Since p e increases in 
comparision to the ideal case, p £ = 2|e|/[(71u) 2 7r], the 
energy-dependent renormalized velocity, vV £ decreases 
up to 10% if g <0.5 and the concentration of carriers 
in an intrinsic graphene increases up to 2 times, see Figs. 
2b and 2c, respectively (here Ut — 8.1 x 10 10 cm~ 2 is 
the equilibrium concentration at room temperature and 
at g ->-0). 

Further, we use the standard relation A £tT> = 
i (Gf p - G^ p +] /2?r and transform the probability (9) 
taking into account the energy conservation law: 



W(e,e')^{\C q \ 2 (N q + l)} r 
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FIG. 3: Contour plots of dimensionless kernels w(e/T, e /T) 

for g =0.4 (a) and g =0.2 (b). 
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The trace here should be taken using tr I PL 

[l+(p-p')/pp']/2andtr(^ ±) J P^ ) ) = [l-(p-p')/pp'}/2. 

This result differs from the standard consideration, |14l | 
because interference of electron and hole states gives an 
essential contribution to W(e, e') due to the matrix struc- 
ture of the spectral density functions. After the integra- 
tions over p-plane, one transforms (11) into 

@GRw(e/T,-e'/T)/T 2 , 



W(e,e') 
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where we separated the dimensionless kernel, w(£, £'), 
and the factor, Qgr, which is written for the case of 
the deformation interaction of carriers with the in-plane 
acoustic modes, see Refs. 10 and 15. Here D is the de- 
formation potential, s is the sound velocity, and p s is the 
sheet density of graphene. At room temperature and typ- 
ical other parameters [lj| we obtain v ac ~ 0.96 x 10 6 cm/s 




FIG. 4: (Color online) Averaged kernel w g versus coupling 
constant g. 

and Ogr — 5.06 x 10 19 cm _2 s _1 (notice, that v ac ex D 2 
and we used D ~ 12 eV). The dimensionless kernel is 
plotted in Fig. 3 and the probability W(e, e') is sup- 
pressed fast if (e — e')/T >0.15. This cut-off factor is de- 
termined by the weak ratio s/v ~1/137 mainly while pa- 
rameter g determines a peak value of W(e,e'), c.f. Figs. 
3a and 3b. 

The generation-recombination rate (6) is written 
through (11) and (12) with the use of the dimensionless 
variables £ = e/T and £' = s'/T: 



„-£'-»> /T 



dU d£'w(£,?) (13) 
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For typical concentrations of carriers, fxf and pf exceed 
0.15T and one can simplify the rate as follows: 
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(14) 



where the averaged over energies kernel w g is plotted ver- 
sus g in Fig. 4, together with a simple parabolic fit. For 
the disorder level corresponding to the resistance ~5 kSl 
per square, one obtains w g ~0.02. 



III. RESULTS 

In this section we analyze the concentration balance 
equation dnt/dt — (dn/dt) ac , where the right-hand side 
of Eq. (14) is written through ipf" = pf jT and ipf~ = 
pf /T, together with the initial conditions ipf =0 = ipQ 
and ip^=o = V'o" determined through the initial concen- 
trations n t= o and nt=o according to Eq. (4). Variables 
ij}f and ipf are connected through the electroneutrality 
condition 



dep e 



1 



3 e/T-4>> + 1 e e/T+i>< + 1 



(15) 



and below we consider the cases of an intrinsic graphene 
(n s = 0) and a n-type heavily-doped graphene (n s > 0). 



A. Intrinsic graphene 



For the case under consideration, ipf = —ipf- = ipt 
and the concentration of electrons (or holes, because now 
n t = n t ) is given by n t = f£° dep e [exp(e/T - ip t ) + 1]~" , 
so that n t a- n( i ipt are connected through 



dn t 
~dt 



dipt 



dep e 



dt J 1 + cosh (e/T - ip t ) ' 
o 



(16) 



As a result, the concentration balance equation is derived 
from Eq. (13) as 



dn t 
dt 



-@GRWg tanh 



</>* 



(17) 



and Eqs. (16) and (17) are transformed into the first- 
order differential equation for ipt with the initial condi- 
tion ipt=o = ipo where ipo is determined through the n{=o- 
The implicit solution of this equation takes form: 



V>o 



vgr% = / dipF(tp), v GR = ir. 
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Here rj = P£t/7>t i s the dimensionless density of states 
and the temporal evolution of n t is described through 
the characteristic rate vgr an d the dimensionless func- 
tion F{ip). At room temperature and at w g ~0.02, one 
obtains vqr — 1-85 x 10 7 s _1 for the parameters used. 

Figure 5 shows the transient evolution of nt normal- 
ized to the equilibrium concentration n eq , for the cases 
of recombination or generation of carriers, if n t= o > n eq 
or n t= o < n eq , respectively. The relaxation becomes sup- 
pressed if the disorder level decreases both due to a slow- 
ness of dependency on VQRt, c.f. curves for g =0.5 and 
0.25 in Fig. 5, and, mainly, due to the relation vqr oc w g , 
see Fig. 4. Within a 5% accuracy, the evolution of n t can 
be fitted by the exponential dependencies 



n t » n + n T [1 - exp(-av GR t)] 



(19) 



with the parameter a varying between 0.16 and 0.28 de- 
pending on the initial conditions. The corresponding 
times, {avQn)~ l , vary between 310 and 180 ns for g ~0.5. 
This time scale is comparable to the radiative recombi- 
nation times, see Appendix. 



B. Heavily-doped graphene 



We turn now to the case of a heavily doped graphene, 
when ipt ^$> 1 and it is convenient to introduce a weak 




FIG. 5: (Color online) Transient evolution of concentration n t 
at different initial conditions: n t= o — 3n eq (1,2), n t= o = 2n 



0.5n eq (5,6) for coupling parameters g 



eq , 

0.25 



(3,4) and n*=o 

(1, 3, 5) and g =0.5 (2, 4, 6). Dotted curves correspond to 
exponential fits (19), with a =0.16 (1), 0.125 (2), 0.19 (3), 
0.155 (4), 0.22 (5), and 0.275 (6). 



variation 5ipt = ipf 



— ip s where ip s corresponds to the 
equilibrium case. Neglecting a hole concentration and 
using the step function in c-band, one obtains n s from 
Eq. (15): 



i>sT 



depe 



(20) 



As a result, 5ipt and ip^ are connected by the electroneu- 
trality condition (15) as follows 



Ht 



Pt 



din 



Pe=^ 3 T J 1 + exp (£ + i>f ) ' 



(21) 



where the ratio ~p T / p e= ^, 3 T can be found from Fig. 2a. 
Using Eq. (14) we transform the concentration balance 
equation into the form: 



d&tpt 
dt 



vgrPt 



2 / o e= , / , sT [l + exp(V> t < )]' 



(22) 



Substituting the relation (21) into Eq. (22) one obtains 
the first-order differential equation for ipf, with the im- 
plicit solution 

VGRt = dtp 



d£rt (1 + eF) 
1 + cosh (£ + ip) ' 



(23) 



i>a 



where ip Q = ip^Lg appears from the initial condition. No- 
tice, that the factor ~p T /p E= ^ sT drops out from the solu- 
tion (23), i. e. the transient process under consideration 
does not depend on the doping level because the only 
low-energy states are involved in the interband transi- 
tions. 




FIG. 6: (Color online) Transient evolution of concentration 
nt — n s at different initial conditions: nt=o — n s = 3tit (1), 
Tlt = o — n s — 2n,T, (2) and nt=o ~ n s — 0.5wt (3). Solid 
and dashed curves correspond to coupling parameters g =0.5 
and g =0.25, respectively. Dotted curves correspond to the 
exponential fits. 



Further, we plot the transient evolution of the hole 
concentration n t = n t — n s determined through ipf ac- 
cording to n t = Jq 00 dep E [exp(e/T + ip<) + 1] . Fig- 
ure 6 shows the concentration n t versus dimensionless 
time, vent, for the initial conditions written through 
tit- Similarly to the undoped case, the exponential fits 
(n t - n s )/(n t=0 - n s ) ss ex.p(-/3v GR t) with (3 ~0.12 (1), 
0.15 (2), and 0.2 (3) describe the transient evolution with 
an accuracy ~10% if vgr^ <10. An enhancement of re- 
combination takes place at tails of transient evolution, if 
v gr% >10. Since the relaxation rate increases with the 
disorder level, vgr °c w g , the recombination process be- 
comes faster in spite of an opposite dependency on vont 
in Fig. 6. The relaxation times, ~ {P^gr) > vary be- 
tween 410 and 240 ns for g =0.5 and different initial 
conditions. Once again, the recombination scale is com- 
parable to the radiative recombination process shown in 
Fig. 7b, Appendix. 



IV. SUMMARY AND CONCLUSIONS 

We have examined the new channel for interband 
generation-recombination process of carriers in a disor- 
dered graphene via acoustic phonons. The efficiency of 
transitions increases with the disorder level and concen- 
tration of nonequilibrium carriers as well as with temper- 
ature. We have found that the relaxation rate belongs to 
submicrosecond range for the samples with typical disor- 
der level at room temperature. 

Let us discuss the assumptions used in the presented 
calculations. The main restriction of the results is the 
description of the response in the framework of the 
quasiequilibrium approach, with different chemical po- 
tentials in c- and w-bands but the same temperature 
due to the fast energy relaxation caused by phonon and 
carrier-carrier scattering processes. We also restrict our- 
selves by the simplest model of the short-range disorder. 



By analogy with the description of transport phenom- 
ena, [l|, [l(J, LL2( more complicated calculations for finite- 
range disorder should give similar results. But the case 
of impurities with a low-energy resonant level, which was 
discussed recently in Refs. 16, requires a special con- 
sideration. We considered the deformation interaction 
of carriers with longitudinal acoustic modes [10 . Il5| ne- 
glecting scattering by surface phonons of the substrate 
in agreement with the experimental data. [17| Such a 
contribution can only restrict the energies under consid- 
eration because of the lower surface phonon energy (^55 
meV for the SiC>2 substrate). Since the non-diagonal 
components give a weak contribution to the concentra- 
tion balance equation under consideration, [111 ] we take 
into account only diagonal components of the density ma- 
trix f at while evaluating of the generation-recombination 
rate. The simplifications mentioned above do not change 
either the peculiarities of the generation-recombination 
processes or the numerical estimates of relaxation times 
given in Sec. III. 

Next, we briefly consider some possibilities for experi- 
mental verification of the mechanism of interband tran- 
sitions suggested. It is clear from a comparison of the 
results in Sec. Ill and in Appendix that interband transi- 
tions via acoustic phonons and via thermal radiation can 
be separated due to different temperature and concentra- 
tion dependencies of damping. A possible contribution 
of the disorder-induced Auger process is beyond of our 
consideration and requires a special study. In contrast to 
the ultrafast optical measurements applied for the study 
of the relaxation and recombination of high-energy car- 
riers, [2J a transient evolution of concentration over time 
scales ~100 ns can be measured directly (e.g. in Ref. 5 
the transient response under abrupt switching on of a dc 
field lasts up to hundreds of nanoseconds) . But under a 
verification of the slow process examined, a possible con- 
tact injection or a trapping into substrate states should 
be analyzed. 

To conclude, we believe that the generation- 
recombination via acoustic phonons can be verified ex- 
perimentally and more detailed numerical calculations 
are necessary in order to separate this mechanism from 
other contributions. The results obtained will stimulate a 
further study of the generation-recombination processes 
which are essential in many transport and optical phe- 
nomena far from equilibrium. 



Appendix: Radiative transitions 

Below we describe the generation-recombination pro- 
cesses which are associated with the interband transitions 
induced by the thermal radiation and evaluate the radia- 
tive relaxation rate for the weak disorder case, j« 1. 
The corresponding collision integral was evaluated in Ref. 
9 and the kinetic equation for the electron distribution 
f e pt takes the form 



9 fept 

dt 



— v p [ N 2vp/T (1 - fept - fhpt) - feptfl 



hpt\ 



(Al) 

where N 2vp /T describes the Planck distribution of the 
thermal photons at temperature T. The hole distri- 
bution can be obtained from the condition d(f ep t + 
fhpt)/dt = 0. The interband absorption or emission 
of photons are described by the first or second terms 
in the right-hand side of Eq. (A.l) and are responsi- 
ble for the generation or recombination processes. The 
rate of spontaneous radiative transitions is given by 
Up — vjip/h where we have introduced the character- 
istic velocity vr ~41.6 cm/s for graphene surrounded 
by SiC>2 layers. Similar to Eq. (6) contribution of the 
radiative collision integral from (A.l) into the concen- 
tration balance equation takes the form (dn/dt) ac — 

(4/£ 2 ) Ep 4 R) [N2vp/T(l - fept ~ fhpt) ~ feptfhpt] ■ 

For the case of an intrinsic graphene, the balance equa- 
tion is written by analogy with Sect. Ill A through 

ipt = I-H/T: 



dipt 
dt 



~VR 



F l {^t) 



VR 



2v R T 



(A.2) 



where v'Z ~30 ns is the radiative recombination time at 
room temperature. The functions Fi(ip) and N(tp) are 
given by 



FiW = 



dtf(l 



,-ty\ 



(l-e- 2 «)(e«^ + ir 



N(j>) = 



dft 



1 + cosh(£ - ip) 



(A.3) 



and the implicit solution of Eqs. (A.2), (A.3) is given by 
the similar to Eq. (18) formula: 



vst 



di/> 



Wo 






(A.4) 
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In Fig. 7a we plot the transient evolution of concen- 
tration versus the dimensionless time, VrI for the same 
initial conditions as in Fig. 5. These transient depen- 
dencies are described by the exponential decay given by 
Eq. (19) with a psO.25 for all cases. Thus, one obtains 
the radiative recombination time ~ A/vr «120 ns which 
does not depend on an initial concentration. 




FIG. 7: (Color online) Transient evolution of concentration 
due to interband radiative transitions for intrinsic (a) and 
heavily-doped (b) graphene at different initial conditions: 
rit=o = 3n T (1), nt=o — 2n T (2) and n t=0 = 0.5n T (3). 
Dotted curves correspond to exponential fits. 



For the case of doped graphene, the concentration bal- 



ance equation (22) should be replaced by 

dStpt 



dt 



fW), 



(A.5) 



FdW 



d^e 



(l_ e -2?)( e e+^ + l) 



while the relation between dipt and ip^ takes form [c.f. 
Eq. (21)] 






d& 



l + exp(S + w) 



(A.6) 



As a result, the equation for ipf nas the only difference 
from Eq. (2) due to the replacement Fi(ip) by Fd(i/j). The 
implicit solution of Eqs. (A.5) and (A.6) is given by (A. 4) 
with the same replacement. In Fig. 7b we plot the tran- 
sient evolution of hole concentration, n t /riT, for the same 
initial conditions as in Fig. 6. The corresponding expo- 
nential fits are determined by the coefficients f$ ~0.15 
(1), 0.125 (2), and 0.06 (3), i.e. the relaxation rate de- 
pends on hole concentration. At room temperature the 
radiative recombination time (/Sz/r)" 1 corresponds to the 
time interval between 190 and 480 ns. 
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